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ABSTRACT 

To characterize the magnetic field structure of the outflow and core region within a prototypical 
high-mass star-forming region, we analyzed polarized CO(3-2) - for the first time observed with the 
Submillimeter Array - as well as 880 /im submm continuum emission from the high-mass outflow/disk 
system IRAS 18089-1732. Both emission features with polarization degrees at a few percent level 
indicate that the magnetic field structure is largely aligned with the outflow /jet orientation from the 
small core scales to the larger outflow scales. Although quantitative estimates are crude, the analysis 
indicates that turbulent energy dominates over magnetic energy. The data also suggest a magnetic 
field strength increase from the lower-density envelope to the higher-density core. 
Subject headings: Stars: formation - Stars: early-type - Techniques: polarization - Techniques: spec- 
troscopic - ISM: jets and outflows - Stars: individual: IRAS 18089-1732 



1. INTRODUCTION 

The influence of magnetic fields on star formation pro- 
cesses has been intensely discussed since decades (e.g., 
ICrutcherl [20051 For magnetically supported molecular 
clouds, the magnetic field structure is expected to de- 
form during the cloud collapse, and theory predicts an 
hour-glass shaped magnet ic field structure pinc hed to- 
ward the protostar (e.g., IStahler fc Pallal 120051) . This 
pinched magnetic field morphology is also required in 
the classical jet-launching scenario, where the gas is 
ejected from the accretion disk along the b e nt magnetic 
field lines (e.g., iBlandford fc Pavnel 119821 : ICamenzindl 
1990). Submillimeter Array (SMA) observations of lin- 
early polarized dust continuum emission toward the low- 
mass protostar NGC1333 IRAS4 have clearly shown 
the pinch ed hour-glass shap ed magnetic field field mor- 
phology (|Girart et all 120061 ). The quest for magnetic 
fields is as important in high-mass star formation. Re- 
cent submm dust polarization observations have revealed 
similarly shaped magnetic field m orphologies toward 
two massive sta r-forming regions (|Girart et all 120091 : 
iTang et "all l2009f) . In addition to that, molecular out- 
flows are known to be as ubiquitous in high- as in low- 
mass star formation, and if the launching mechanisms 
are of similar nature, magnetic fields have to be present 
in massive cores as well. 

While most previous high-spatial-resolution magnetic 
field studies were either based on dust cont inuum polar- 
ization observations or maser studies fe.g.. iGirart et all 
l2009trvTemming s 2008b), the molecular gas can also show 
non-maser spectral line pola rization signatures due to 
the G oldreich-Kylafis effect (jGoldreich fc Kvlafisl Il98ll 
1982). Observations of this effect are particularly promis- 
ing for studying the magnetic field in molecular outflows 
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(IGirart et all 119991 : lLai et all 120031 : ICortes et all I2005L 
120081 ). 

We present a combined polarization study of the 
880 /im submm continuum and the CO (3-2) emission 
toward the well-studied massive disk-outflow system 
IRAS 18089-1732. The luminosity and gas mass of the 
region at a di stance of 3.6 kpc are lO 4 5 L w and 1200 Mq , 
respectively (ISridharan et all 120021: IBeuther et"all l2002t 
IWilliams et all 12004 . IVlemmingsl (|2008al) derives for 
this region a line-of-sight magnetic field strength of 
8.4 mG from Zeeman splitting measurements of the 
6.7 GHz CH3OH maser line at densities > 10 6 cm~ 3 . 
SMA studies revealed a molecular outflow in approxi- 
mately north-south direction and rotational signatures 
in the dense gas perpendicular to the outflow which 
were confirmed by high-excitation NH3 obs ervations 
(|Beuther et all 120051 : IBeuther fc Walshl I2008D . Since 
IRAS 18089-1732 remains a single compact source up 
to the highest angular resolution (< 1"), and since the 
outflow-disk orientation is well constrained, this region is 
the ideal target for magnetic field investigation in high- 
mass star formation. 

2. OBSERVATIONS 

The region was observed with the SMA 5 ) in the com- 
pact and extended configurations in May 2008 and Au- 
gust 2009. The phase center was R.A. 18:11:51.4 and 
Dec: -17:31:28.5 (J2000.0) and the tuning frequency 
345.796GHz in the upper sideband (z;i sr = 33.8 km s -1 ). 
While the compact array observations were still con- 
ducted in the old setup with 2 GHz bandwidth in both 
sidebands, the extended configuration data profited from 
the increased bandwidth of 4 GHz in both sidebands. 
The spectral resolution in both cases was 0.8125 kHz 
(~ 0.7km/s). The weather was good with zenith opac- 
ities mostly below 0.25 at 346 GHz. Calibration an d 
imaging were conducted in MIRIAD (jSault et allll995T ). 
The passband calibration was derived from 3c454.3, and 
fluxes were calibrated via the regularly monitored SMA 
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calibrator database. The flux density scale is estimated 
to be accurate within 15%. For the phase and ampli- 
tude calibration regularly interleaved observations of the 
quasar 1733-130 were conducted. The p olarimeter sys- 
tem o f the SMA is desc r ibed i n detail in iMarrone et al.l 
(2006|); M arrone fc Raol (|2008l ). As polarization calibra- 
tor we used for the compact and extended configura- 
tion 3c273 and 3c454.3, respectively. For the CO(3- 
2) polarization calibration we took only the emission 
from the 7 chunks of the bandpass around the CO (3- 
2) line (approximately 700 MHz). For the continuum 
data we used both configurations resulting in a syn- 
thesized beam of 1.65" x 1.05" (P.A. 51°). Since the 
CO (3-2) emission is more extended, we used only data 
from the compact configuration resulting in a synthe- 
sized beam of 2.11" x 1.34" (P.A. 52°). The la rms val- 
ues of the continuum Stokes I and Stokes Q/U images 
are 18.4 mJy beam -1 and 4.3 mJy beam -1 , respectively, 
because the noise is dominated by side-lobes from strong 
sources due to insufficient uv-sampling. For the CO (3-2) 
emission the la values at 4kms -1 spectral resolution of 
the Stokes I and Stokes Q/U images are 140 mJy beam -1 
and 75mJybeam -1 , respectively. The polarized contin- 
uum and CO(3-2) images were produced applying a 3a 
cutoff to the data. For comparison, for the continuum 
we also show the results applying 2a cutoffs in Fig. [TJ 
From the other spectral lines in the setup, here we only 
present the integrated SiO(8-7) emission (Fig.Q]) to out- 
line the direction of the outflow/jet. The synthesized 
beam and the rms of the SiO(8-7) image integrated 
from 30 to 40kms -1 are 2.11" x 1.33" (P.A. 52°) and 
52 m Jy beam - 1 , respectively. 

3. RESULTS 
3.1. Submm dust continuum polarization 

Figure Q] presents the linearly polarized ( \J Q 2 + U 2 ) 
as well as Stokes I 880 /im continuum emission to- 
ward the outflow/disk system IRAS 18089-1732. While 
the Stokes I image exhibits the same single-peaked 
structure as known from previous observations (e.g., 
iBeuther et al.l I2005D . the linearly polarized emission 
has sub-structure offset from the main peak. The 
strongest linearly polarized emission with a peak flux 
of 23.9 mJy beam -1 is ~ 0.7" N-NE (north-northeast) of 
the Stokes I peak, and the polarized emission shows addi- 
tional substructure south of the Stokes I peak position. 
At the linearly polarized peak position (R.A. (J2000.0) 
18:11:51.51, Dec. (J2000.0) -17:31:28.5), the polarization 
fraction is at a ~1.5% level, at the Stokes I peak position 
(R.A. (J2000.0) 18:11:51.47, Dec. (J2000.0) -17:31:28.8), 
the polarization fraction is still at a ~1.1% level (Ta- 
bled]). The Stokes I peak and integrated fluxes within 
the 10<7 contours are 1.93 mJy beam -1 and 4.24 Jy, at 
880 /ur n wavelength, respectively. Following Hildcbrand 
(1983), assuming optically thin dust continuum emission 
with a dust e missivity index = 2 at a temperature of 
~100K (e.g. . IBeuther et alJl2004l) and a gas-to-dust ra- 
tio of 186 (e.g.. iDraine et al1l2007ft . the derived gas col- 
umn densities N and masses M are 9.5 x 10 24 cm -2 and 
197 M©, respectively. Taking into account that in partic- 
ular the mass is a lower limit due to the interferometric 
missing flux, clearly we are dealing with a high column 
density and massive core. Assuming a spherical distri- 
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Fig. 1. — The grey-scale presents the linearly polarized 880 /im 
continuum image of IRAS 18089-1732 in units of Jy beam - 1 . The 
line-segments plotted in approximate half-beam spacing of 0.5" 
show the magnetic field orientation assuming that it is oriented 
perpendicular to the polarization (line-segment length scales with 
strength of polarized emission). Line-segments marked by a small 
green circle are those where a 2a cutoff was applied. The red con- 
tours present the integrated Stokes I image of the continuum emis- 
sion in lOcr steps. To outline the outflow direction, the blue con- 
tours show the SiO(8-7) emission integrated from 30 to 40kms — 1 
(5cr steps of 260 mJy beam - 1 ). Full and open ellipses present the 
synthesized beams of the continuum and SiO emission. 

bution with an approximate core diameter of 3.5" (linear 
size of ~ 12600 AU), we get an approximate average core 
density n of 5 x 10 7 cm -3 . 

For elliptical dust grains, the polarization orientation 
of the radiation field is expected to be perpendicular to 
the magnetic field if the dust grains spin fast enough. 
Therefore, to infer information about the magnetic field 
structure, we rotated the line-segment orientation of the 
polarized emission in Figure [1] by 90° (line-segments are 
shown in approximate Nyquist sampling of 0.5", indicat- 
ing the number of independent measurements). In this 
picture, the overall orientation of the magnetic field com- 
ponent in the plane of the sky follows approximately the 
general morphological structure of the polarized emis- 
sion with an orientation roughly along an axis in N-NE 
to S-SW direction. This general orientation is approxi- 
mately aligned with t he outflow traced by the SiO (8-7) 
emission (Fig. Q] and IBeuther et al]l2004h . Toward the 
western edge of the core, the field orientation bends a 
little bit toward the mid-plane of the underlying accre- 
tion disk HMt^&^iSllQQl). 

3.2. Polarization of the CO (3-2) outflow emission 

Furthermore, we observed simultaneously the polar- 
ized emission of CO(3-2). Figure [2] presents the Stokes 
I and the linearly polarized Stokes U and Q spectra ob- 
tained vector-averaged on the shortest baseline. While 
the emission at the vi sr is filtered out in all cases, the ve- 
locity component around 40kms -1 is not only detected 
in Stokes I but also in both linearly polarized compo- 
nents of Stokes U and Q. From these data, we derive spa- 
tially and spectrally resolved images of the total linearly 
polarized CO (3-2) emission in this region (\/Q 2 + U 2 ). 
Figure [3] shows two 4kms -1 wide channels of polarized 
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Fig. 2. — The three panels present vector-averaged spectra of 
the CO(3-2) line taken on the shortest baseline of the Stokes I, Q 
and U components. The vertical line marks the Uj sr . 

CO(3-2) emission around the velocity of 40kms _1 . We 
again clearly identify the overall N-NE to S-SW outflow 
direction, but here not only the northern component but 
also a southern counterpart. Since we do not identify 
a clear blue-red-shifted dichotomy, the outflow is likely 
close to the plane of the sky. The orientation of the po- 
larization follows mainly the outflow direction, and we 
find polarization degrees in the images of the order of 
10% (Table [J), higher t han found toward other regions 
(e.g., ICortes et al 
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3.3. Magnetic field strength 

Although magnetic field estimates based on polariza- 
tion data are relatively unreliable since they are based on 
a statistical analysis of the polarization line-segments, 
nevertheless, we have enough independent polarization 
measurements in the continuum and line emission that 
at least rough estimates are useful. 

Dust continuum polarizati o n — F ollowing the work 
by IChandrasekhar fc Fermil (|1953l ) and more recently 
iGirart et al.l (|2006l ). we estimate the strength of the 
magnetic field component in the plane of the sky 



B 



with p 
Av los 



pos ~ ^v^"P A0 

rived above in gcm~ 3 
locity dispersion along the line of sight (Av\, 
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the average density de- 
the approximate ID ve- 



2.2 km s" 1 with Au(CH 3 CN) 



At;(CH 3 CNWin2_ 

5.1 km s -1 , iBeuther et al.l l2005i) and A<j> the approxi- 
mate flux-weighted dispersion of the polarization angle 
line-segments of ^27°. In practice, we use the rms 
of flux-weighted line-segments ignoring any bending of 
the field lines which implies that the inferred B-values 
should be lower limits. Q is a dimensionless param- 



eter accountin g for the cloud structure, and work by 
lOstriker et al.l (|2001l ) suggests Q ~ 0.5 for turbulent 
clouds. Since the density p depends on the distance 
d, and the latter is the relatively uncertain kinematic 
distance, one should keep in mind an additional in- 
dependence of B oc 1/ -\/c?/3.6kpc. The derived magnetic 
field str ength in the plane o f the sky is B pos w 11 mG. 
Because IVlemmingsl ([2008a) derives a comparable line- 
of-sight magnetic field strength for similar densities of 
the gas from maser Zeeman observations (B\ os ~ 8mG), 
the total magnetic field strength at the given core density 

should be B to t f 



B 2 

pos 



B \os 



14 mG. 



While this is larger than the value previously derived 
in DR21 (1 mG, albeit meas ured at a density more than 
an order of magnitude lower. lLai et alj|2003t) . it is of the 
same order as the recent estimates based on dust po- 
larization measurements for G31.41 (jGirart et al.l I2009T) 
and those derived from CH 3 OH maser Zeeman measure- 
ments to wards a sample of h igh-mass star- forming re - 
gions by IVlemmingsl (|2008al ); IVlemmings et all (|2010t) ; 
iSurcis et al.1 (|2009f ). The corresponding Alfvenic velocity 
is va oc B I \J~p~ « 3.1kms _1 . The important mass-to-flux 

ratio can be estimated via M/$ B ~ 7.6 x 10~ 24 ^§^ 
5.2 in units of the critical mass-to-flux ratio ( M/<l>B) cr it 
(iVfHg) and B are given in cm " 1 and mG, ICrutcherl 
119991 : iTroland fc Crutcherl 120081 ) . While the latter es- 
timate indicates that the magnetic field cannot pre- 
vent the core from collapse, comparison of the Alfvenic 
velocity with the ID velocity dispersion suggests that 
the collapse proceeds at approximately or s lightly below 
the A lfvenic speed. Furthermore following IGirart et al.l 
we estimate the ratio of turbulent to magnetic 
energy via j3 « 3(Av\ os /va 1ob ) 2 ~ 2.5 (with va 1o , ~ 
2.4kms _1 the Alfvenic velocity along the line of sight). 
Because the estimated B pos is a lower limit, the ra- 
tios M/$b and /3 are upper limits. While for G31.4, 
W75N, W51 and CepA the magnetic field may play a 
dominant role determining the dynamics of the systems 
(IGirart et aI]|20q|jSurcis et alj 120091 iTang et all 12001 
IVlemmings et al.ll2010D . for IRAS 18089-1732 the system 
appears to be more strongly dominated by turbulent mo- 
tions. 

CO(3-2) polarization — Following the same 
Chandrasekhar-Fermi approach as above, we derive 
a very rough estimate of the magnetic field strength 
within the outflow lobes. The rms of the polarization an- 
gle distribution in the northern lobe is A<f) ~ 45° (again 
measured as the rms of flux-weighted line-segments 
ignoring any bending of the field lines). However, 
getting a density estimate based on the given data is 
far more difficult. Since we cannot derive the density 
ourselves, we can only assume an average density in the 
more diffuse envelope where the outflow is observed. Av- 
erage densities of order 10 3 cm -3 should approximately 
represent the envelope/outflow. With these values and 
the same line width as for the core, we get magnetic 
field estimates of order ~ 28 pG, about three orders of 
magnitude below the central core magnetic field values. 
While the Chandrasekhar-Fermi method in principle 
is prone to many uncertainties - e.g., the A<f> requires 
subtraction of the ordered magnetic field - it may even 
be less applicable in the case of the outflow lobes since 
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Fig. 3. — The grey-scale presents the linearly polarized CO(3-2) emission in two 4kms— 1 wide channels in units of Jybeam - 1 (central 
velocities are marked in the top- left corners). The line-segments plotted in approximate half-beam spacing of 0.75" show the orientation 
of the polarized emission (line-segment length scales with strength). The contours present the corresponding CO (3— 2) Stokes I image in 
lOtr steps (positive and negative features as full and dashed contours). The cross marks the position of the submm continuum peak. The 
synthesized beam is shown at the bottom left of each panel. 

the outflow may determine most of the polarization 
orientation. Furthermore, the density uncertainties add 
to the total uncertainty of the magnetic field estimate. 
While we cannot give a reliable estimate of this total 
uncertainty, it is nevertheless likely below 3 orders of 
magnitude. Hence, the magnetic field strength in the 
outflow region appears lower than in the core region. 
This is consistent with previous results for DR21, where 
the CO polarization traces magnetic fields of the order 
10 fiG in the envelope and the dust polarization field 
strengths of the o r der mG in the denser core regions 
(jCortes et al.ll2005HLai et al.ll2003t ). 



4. DISCUSSION AND CONCLUSION 

The combined polarization data of the submm con- 
tinuum and the CO(3-2) emission allow us to study 
the magnetic field structure in the core and outflow 
of the high-mass disk-outflow system IRAS 18089-1732. 
For the dust continuum emission, compared to mag- 
netic field hour glass sha pe structures like t hose found 
in NGC1333 or G31.41 (jOirart et all l2(M I2009T ). the 
case is less clear here. Inferring an hour-glass shape like 
structure indicative of ambipolar diffusion from these 
data would likely be an over-interpretation. However, 
in a br oader sense this may not b e suc h a surprise. 
Already iMouschovias fc Paleologoul (|1979f) pointed out 
that adding rotation to the collapsing core could destroy 
the hour-glass shape structure. Therefore, in a statisti- 
cal sense it may be that finding the hour-glass signatures 
could be more the exception than the rule. An analysis of 
a larger sample of sources is required to confirm or reject 
this hypothesis. Nevertheless, the data of IRAS 18089- 
1732 are consistent with a picture where the outflow and 



magnetic field orientation are tightly linked. 

For the CO (3-2) emission, since the Goldreich-Kylafis 
effect occurs only when the line emission is dominated by 
radiative excitation and not collisions, it sensitively de- 
pends on the anisotropy of the optical depth and the ra- 
diation field as well as on the orientation of the magneti c 
field (for a qualitative descriptions see iKvlafisI 119831) . 
Additionally, there exists an ambiguity whether the po- 
larization orientation is perpendicular or pa rallel to the 
magnetic field. However. ICortes et al.l (|2005j ) have shown 
that a large radiation and optical depth anisotropy pro- 
duces high degrees of polarization parallel to the mag- 
netic field (their Fig. 7). Since we observe such a high 
polarization degree of order 10%, and since the radiation 
and velocity field is dominated by an anisotropic outflow 
cone with a central protostellar radiation source, this sce- 
nario supports an interpretation that the magnetic field 
and CO(3-2) polarization orientation in IRAS 18089- 
1732 are spatially aligned. This is consistent with pre- 
vious CO polariza t ion studies in D R21 and NGC1333 
(jGirart et al.lll999t lLai et al.l I2003T ) as well as with the 
picture of magnetic field accelerated and collimated jets 
where the magnetic field and outflow should have the 
same orientation. 

To the authors' knowledge, this is the first successful 
detection of the Goldreich-Kylafis effect in the CO (3-2) 
line. Comparing the outflow morphology with the po- 
larization direction, it is likely that the magnetic field 
and polarization both are aligned with the outflow. Al- 
though the Chandrasekhar-Fermi method allows only 
rough magnetic field strength estimates, the data are in- 
dicative of a magnetic field strength increase from the 
envelope toward the core. Furthermore, the observations 
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TABLE 1 

Additional online table: Polarization measurements with respect to 0/0 position R.A. (J2000.0) 18:11:51.40 Dec. (J2000.0) 
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